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ABSTRACT 



We derive a simple closed- form expression, relating v 12 (r) - the mean relative velocity 



of pairs of galaxies at fixed separation r - to the two-point correlation function of mass 
density fluctuations, We compare our analytic model for v 12 (r) with N-body 

q ' simulations, and find excellent agreement in the entire dynamical range probed by the 

simulations (0.1 ~ £ ~ 1000). Our results can be used to estimate the cosmological 
density parameter, 0, directly from redshift-distance surveys, like Mark III. 



Subject headings: Cosmology: theory - observation - peculiar velocities: large scale 
J> ' flows 
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■ 1. An analytical model for v 1 

OO ' 

Most dynamical estimates of the cosmological density parameter, f2, use the gravitational effect 
• of departures from a strictly homogeneous distribution of mass. One such dynamical estimator 

can be constructed by using an equation expressing the conservation of particle pairs in a self- 



gravitating gas. This equation was derived by Davis & Peebles (1977)1 from the BBGKY theory 



(see also [Peebles 1980| , hereafter LSS). Consider a pair of particles at a comoving separation vector 
x and cosmological time i, moving with a mean (pair-weighted) relative velocity v 12 (x,t)x jx. It's 
magnitude is related to the two-point correlation function of density fluctuations, £(x,i) , by the 



pair conservation equation (LSS), 



c3 1 a d£{x,a) _ _ v 12 (x,a) , , 

3[l+£(x,a)] da ~ Hr ' U 

where a(t) is the expansion factor, r = ax is the proper separation, H (a) is the Hubble parameter, 
while £(x,a) is the two-point correlation function, averaged over a ball of comoving radius x : 

f(x,a) = 3x~ 3 f Z(y,a)y 2 dy . (2) 
Jo 

At the present cosmological time a = 1, x = r and H = 100 /ikm s _1 Mpc~ 1 , where h is the 
conventional dimensionless parameterization for the Hubble constant. There are two well known 
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approximate solutions of ([!]). They are: the small separation limit, where 3> 1 (stable clus- 
tering regime), and the large separation limit, where |£| <C 1 (linear regime). The stable clustering 
solution is (LSS, §71) 

v 12 (x,a) = -Hr , (3) 

as expected for virialized clusters on sufficiently small scales. The linear solution is given by the 
first terms in perturbative expansions for v 12 and £, which for the correlation is given by 

c = + e } +... , (4) 

with £ (2) = O [£ (1) ] 2 , etc. The growing mode of the linear solution is 

£«(x,a) = Ux)D(a) 2 , (5) 

where D{a) is the usual linear growth factor (LSS, §11; we neglect the decaying mode). The general 
technique for deriving £ (2) and higher order terms for initially Gaussian density fluctuations in an 
Einstein-de Sitter universe was introduced by Juszkiewicz et al. (1980, 1984), Vishniac (1983) and 
Fry (1984); their results were recently generalized to a wider class of cosmological models including 
those considered below ( Juszkiewicz et al. 1993 , Bouchet et al. 199E| ). These calculations show 



that the second order term in the expansion for £ can be written as 

e\x,a) = &x)D(a)*, (6) 

where £ 2 is a function of x alone. Substituting £ = D 2 £ x + D A £ 2 into eq. (||) and solving for v 12 

we get 

v 12 = -iHrf [£W - + 2^] + O^f , (7) 

where / = din D/ din a ~ Q ' 6 (which is a good approximation if A = or Q + Q\ = 1; S7a = 
A/3H 2 ). The bars above £-s denote averaging over a ball of comoving radius x. If the logarithmic 
slope of the correlation function, ~f(x) = — d In £ (1) (x, a) jd In x, is a constant in the range < 7 < 2. 
Then £ (2) is related to by a simple closed-form expression, 

e\x,a) = a( 7 ) [P(x,a)] 2 , (8) 

where a is a constant ( |Lokas et al. 19*9^ ), and the function 0(7) can be expressed in terms of 



Euler's T functions (Scoccimarro & Frieman 1996). In the range < 7 < 2, this expression is well 



approximated by a fitting formula 

a = 1.84 - 1.1 7 - 0.84 ( 7 /2) 10 . (9) 

For 7 > 2, perturbation theory diverges, but N-body simulations suggest that the relation (|8|) 
remains valid, provided the 0(7) dependence for 7 ~ 1.5 is derived from the so-called universal 
scaling relation (USH) - an empirical formula for the nonlinear power spectrum, extracted from 
7 > 1 scale- free and cold dark matter (CDM) simulations ( Jain et al. 1995| , |Scoccimarro fc Frieman 



1996). Since the USH approach fails when 7 < 1, while the perturbation theory (PT) fails for for 



7 > 2, we propose to trade accuracy for an extended range of validity and interpolate the expression 
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ft (7) between the PT formula for 7 < 1 and the USH result for 7 > 1 ( JScoccimarro fc Frieman 
|1996| , Fig. 20), and replace eq. (0) with 



1.2 - 0.657 • 



(10) 



The allowed range of slopes for the above expression is < 7 < 3. To extend the range of validity 
of our perturbative solution (|7|) into the nonlinear regime, we propose the following Ansatz: 



v 12 (x,a) 



I Hrf£(x,a) 1 + ft£(x,a) 



£(x,a) = £(x,a) [1 + £(x,a)] . 

The expression (|ll]) agrees with the perturbative expansion (^) when £ 
the stable clustering limit when x — > 0. Indeed, when £ is large, —v 12 /Hr 
which is of order unity for all reasonable values of f2 and a. 



(11) 
(12) 

0; it also approximates 
(2/3)/(0)(l + «) 



Eq. (11) can be used to predict the relative velocity v 12 at any separation r. The only input 
necessary for that is the correlation function in an interval of separations < r. This is very different 
from the semi-analytic expression for v 12 (r), derived by Mo at al. (1997). In their case, the 
knowledge of the present £(r) is not sufficient to calculate v 12 (r); they use the USH approach and 
a large array of parameters instead. 



2. N-body simulations 

In this section we compare our analytic expression for v 12 with the results from high-resolution 
AP 3 M simulations of 256 3 dark matter particles, kindly provided to us by the Virgo collaboration 
( Jenkins et al. 199§| ). We consider four members of the CDM family: an open model (OCDM), 



a zero curvature low-f2 model (ACDM), and two models with an Einstein-de Sitter metric - the 
'standard CDM' model (SCDM) and its modified version (rCDM). Following Jenkins et al., the 
values, assigned to parameters (h, O, Q A , cr 8 ) are: (0.7, 0.3, 0, 0.85) for OCDM, (0.7, 0.3, 0.7, 0.9) 
for ACDM, and (0.5, 1, 0, 0.6) for the SCDM and rCDM which has extra large-scale power (added 
in an ad hoc manner, described by Jenkins et al. 1998] ) . Here and below a 8 is the rms dark matter 



density contrast in a ball of radius 8 h x Mpc. 

Since CDM-like models are not scale-free, eq. (|H]) does not apply. In principle, we should 
therefore calculate ot(x) = £ (2) (x, a)/[£ (1) (x, a)] 2 for each considered power spectrum and each sepa- 
ration x, using standard perturbative techniques (|Lokas et al. 199*6 , jScoccimarro fc Frieman 1996 ). 



However, as we will show below, these calculations can be significantly simplified by finding an 
effective slope, 7 e ff, which provides 'best fit' a and v 12 when substituted in eq. (10) and (|TT|). The 



precise value of a is unimportant in the stable clustering regime as well as in the linear regime, when 
the term quadratic in£ is sub-dominant. Hence, the precision in a(x) matters only at the boundary 
between the linear and nonlinear regimes, say at £(x,t) = 1. One of the possible definitions of the 
effective slope is therefore given by 

7cff = -(dhi£/dlnz)|£ = i . (13) 
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One can also choose 7 e g as follows. In Fig. 1, we plot logarithmic slopes of and £(x) (the 

latter is measured from the simulations). Both curves agree at large separations as they should, 
apart from small differences arising from noise in the measurement (we use only a finite number 
of bins and pairs to measure £), and from finite box-size and cosmic- variance effects flJenkins~et 



al. 1998| ). However, there is a well-defined scale (marked with an asterisk) at which the non-linear 



slope turns away from the linear theory prediction, marking the onset of the non-linear regime. 
We take 7 c g to be the logarithmic slope of at that scale. The resulting slopes are: 7 e ff = 1.67 
(SCDM), 1.46 (ACDM), 1.49 (OCDM), and 1.40 (rCDM). The alternative definition (|T|) gives, 
respectively: 7 e g- = 1.67, 1.47, 1.45, and 1.28. The advantage of the former definition is that it 
is more closely related to observations because it uses £,(x) only rather than £ along with £ (1) (x). 
The linear correlation function is easy to calculate under the controlled conditions of an N-body 
experiment but it can not be easily determined from observations. 

In Figure 2 we test our Ansatz (|TT|) against N-body measurements. For comparison, we also 
plot three other approximations for v 12 (r), considered earlier in the literature: 

A) v 12 = -lHrft m , B) v 12 = -\Hrfl , C) v 12 = -§#r/f . (14) 

(A) and (B) are two variants of linear theory predictions, and (C) is an improvement over linear 
theory, suggested by Peebles (LSS, §71). Figure 2 shows that the deviations from linear theory 
are small at large separations, as they should. The range of validity of the Peebles formula (C) is 
already considerably wider than that of linear theory. However, our new nsatz provides by far the 
best approximation. In fact, it covers the entire dynamical range! 

The scale at which the linear approximation becomes acceptable depends on the amplitude of 
fluctuations; it increases with increasing cr 8 . For example, for a power-law correlation function, we 
have 

£(r) = cr 8 2 ^(7) r~ 7 , where (15) 
F(i) = (16/i" 1 Mpc) 7 (4-7)(6-7)/24 , (16) 



and eq. (14 B) can be rewritten as 



v 12 {r) « -|(T 8 2 O a6 i/F(7)r 1 - 7 « - 667 a 8 2 ft a6 km/s , (17) 

where the expression after the last "~" sign assumes 7 = 1.75 and r = 10/i -1 Mpc. The relative 
error in the latter expression, introduced by linear theory, can be calculated from eqs. ( |i"T|) and 
(|i~0|). It depends on a$ only; for <r 8 = 1 and 0.6, linear theory overestimates |u 12 | by 24% and 10%, 
respectively. 



3. Velocity bias 

So far we considered the dynamics of pairwise motions of dark matter particles. However, for 
practical applications, it is necessary to understand the relation between v 12 (r) and the relative 
pairwise velocity of galaxies, v 12g (r). We define the galaxy clustering bias as the square root of the 
ratio between the galaxy and the dark matter correlation functions: b(r,t) 2 = £ g (r, i)/£(r, t). 
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Observations suggest that there is no velocity bias: splitting the Mark III catalogue (Willick 



et al. 1997[) into subsamples of elliptical and spiral galaxies has no effect on v 12g flJuszkiewicz et 



al. 1998j )7Similar results follow from recent numerical simulations, which account for dissipative 
processes, important for galaxy formation. These models do show some clustering bias on small 
scales; however, there is no velocity bias, and v 12 = v 12g ( |Kauffmann et al. 1998| ). 

As we show now, this contradicts the simplest toy model for bias, where b ~ l/cr 8 is a constant. 
In this prescription, also known as linear bias, the galaxy density contrast at position r A is simply 
given by 5 gA ~ bS A , where 5 A = Pa/{p) — 1 is the mass density fluctuation, and A = 1,2,... 
enumerate galaxy positions. In the fluid approach, the mean pairwise velocity between two points 
at distance r is given by 

Vl2[r) - TTW) ' [ ' 

where v A is the peculiar velocity at a point f A , r = \r x — r 2 \ is the separation, and £(r) = (J^). 
For the galaxy pair density- weighted relative velocity, v 12g , the matter density field in the above 
expression, <5, has to be replaced by 5 g . In the limit of large separations (5, £ — > 0), the linear 
biasing model, applied to eq. (|lS|), gives v 12g (r) = bv 12 (r), and since v 12 oc one obtains 

v 12g oc cr g ft a6 (Fisher et al. 1994). On small scales, where 1 + S ~ 5 and 1 + £ ~ £, the factors of 
b in the denominator and numerator cancel and Vi 2g (r) = v 12 (r). This unphysical behavior shows 
the limitations of the linear biasing model. Similar results can be obtained from considerations, 
based on the continuity equation and gravitational clustering. Gravity tends to remove bias and 
b has to evolve with time and separation ( |Fry 199~6| ). The results of Kauffmann et al. (1998) 
fit neatly into this picture: their measurements of b(r) at the 'present time' (redshift = 0) are 
significantly different from unity only on small scales where £ ^> 1 and gravity is no longer the only 
force determining the dynamics; at large scales, however, gravity dominates and 6 = 1. Since the 
dependence on b in eq. (|i~8|) cancels out when £ ^ 1, it is not surprising that the simulations give 
Vi 2 = v 12g on all scales. Summarizing, we conclude that a realistic model of biasing should lead 
to unbiased pairwise velocities; the linear bias model disagrees with observations as well as with 
simple theoretical considerations. 



Conclusions 



The main result of this paper is the expression ( pT[ ) - an approximate analytic solution of the 
pair conservation equation. Our equations (|K]) - ( Jl3| ) reproduce the results of numerical simulations 
on all scales from the strongly non-linear to the linear regime. 



Our results can be used to estimate ft from redshift-distance surveys, like Mark III ( [Ferreira 
et al. 1998|, |Juszkiewicz et al. 1998D 



On very large scales, v 12 (r) is proportional to ft°' 6 <7 8 2 . 



On 



intermediate scales (the mildly non-linear regime) this degeneracy is removed and ft and <r 8 can be 
measured separately because the a£ term in eq. ( |IT|) is ft-independent. These properties make our 
estimator complementary to other estimators. Indeed, the POTENT method ( |5igad et al. 1998 ) 
and the cluster abundances ( [Bahcall fc Fan, 199S , |Eke et al. 1998 ) are sensitive to (3 = ft°' 6 cr 8 ; the 
supernovae ( Riess et al. 1998 , perlmutter et al. 1998|) distances measure q Q = (ft/2) — ftA; finally, 
the position of acoustic peaks in the CMB power spectrum ( Doroshkevich et al. 1978[ ) is sensitive 
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to O + J7a- The advantage of our estimator over the CMB peaks method is model-independence; 
its advantage over POTENT is simplicity ( [Ferreira et al. 199Sj ). 

Our results can also be used to study the nature of biasing and to test the gravitational 
instability theory. Indeed, one can use the galaxy correlation function, estimated from galaxy 

redshift or angular surveys to predict v 12 (r) which could then be compared to a v 12 (r) estimated 
from redshift-distance surveys flGaztahaga fc Juszkiewicz 1998 ). 

We thank Antonaldo Diaferio, Enrique Gaztahaga and Roman Scoccimarro for insightful 
comments. This work was supported by grants from the Polish Government (KBN grants No. 
2.P03D.008.13 and 2.P03D.004.13), the Tomalla Foundation, by the Poland-US M. Sklodowska- 
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r [ h' 1 Mpc ] r [ h' 1 Mpc ] 

Fig. 1. — Logarithmic slope din £ / din r of the linear theory correlation function (dashed), and 
the measured non-linear correlation function (solid) for the four Virgo simulations which we have 
analyzed. The asterisks mark the effective slopes 7 e fr used in equations (pi]) and (|i~l]), respectively. 
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Fig. 2. — The mean pairwise velocity v±2 of the Virgo simulations (solid lines) compared with four 
closed-form approximations to solutions of the pair conservation equation: the two versions of the 
linear approximation, eq. [14] (A), and (B), are plotted as thin and thick dashed curves respectively; 
the Peebles approximation, eq. [14| (C), is shown as dot-dashed curve; and eq. ( |TT| ) - the Ansatz 
proposed in this paper - is drawn as dot-dot-dot-dashed line. 



